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Abstract 

Assuming the definition of the inversion problem (IP) as the exact matching of the terms in the low 
redshift expansion of cosmological observables calculated for different cosmological models, we solve 
the IP for Dl(z) and the redshift spherical shell mass density mn{z) for a central observer in a LTB 
space without cosmological constant and a generic AC DM model. We show that the solution of the 
IP is unique, corresponds to a matter density profile which is not smooth at the center and that the 
same conclusions can be reached expanding self-consistently to any order all the relevant quantities. 

On contrary to the case of a single observable inversion problem, it is impossible to solve the IP 
(LTB vs. ACDM) for both mn(z) and Dl(z) while setting one the two functions k(r) or t b (r) to 
zero, even allowing not smooth matter profiles. Our conclusions are general, since they are exclusively 
based on comparing directly physical observables in redshift space, and don't depend on any special 
ansatz or restriction for the functions defining a LTB model. 
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I. INTRODUCTION 

High redshift luminosity distance measurements jl-6] and the WMAP measurements [7, 8| 
of cosmic microwave background (CMB) interpreted in the context of standard FLRW cosmo- 
logical models have strongly disfavored a matter dominated universe, and strongly supported 
a dominant dark energy component, giving rise to a positive cosmological acceleration. As 
an alternative to dark energy, it has been proposed {9, 10] that we may be at the center of 



an inhomogeneous isotropic universe without cosmological constant described by a Lemaitre- 
Tolman-Bondi (LTB) solution of Einstein's field equations, where spatial averaging over one 
expanding and one contracting region is producing a positive averaged acceleration a^, but 
it has been shown how spatial averaging can give rise to averaged quantities which are not 



observable 



Hi 



Another more general approach to map luminosity distance as a function of redshift Dl(z) 
to LTB models has been recently proposed {12I [lsI ] , showing that an inversion method can be 
applied successfully to reproduce the observed Dl(z). 

Interesting analysis of observational data in inhomogeneous models without dark energy 



and of other theoretically related problems is given for example in Our definition of 

the inversion problem is purely mathematical, and consists in matching exactly the terms of 
low redshift expansion for the relevant cosmological observables for different models. Following 
this procedure we determine the LTB models which are mimicking at low redshift the ACDM 
theoretical predictions. 

It is important to observe that according to our definition of the IP we look for the LTB 
models which match the coefficients of the redshift expansion of the observables corresponding 
to the best fit flat ACDM models, but this kind of approach is not completely rigorous since 
while these models depend only on the two independent parameters Hq, Qm, LTB models have 
a higher number of parameters, implying that they could actually provide an even better fit of 
experimental data. From an observational point of view it would in fact be more important to 
confront directly with actual experimental data rather than reproduce the best fit theoretical 
ACDM model, so our conclusions should be considered keeping this limitation in mind, and 
we leave to a future work such a direct analysis of experimental data. Our definition of the IP 
appears nevertheless the most natural one from a mathematical point of view since it consists 
in matching the theoretical predictions of different models within the range of validity of the 
Taylor expansion, and the formulae derived could actually be used even for the purpose of low- 
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redshift data fitting, since they have the advantage of not being dependent on any functional 
ansatz. We could for example obtain contour plots for the coefficients of the expansion of the 
functions k(r) , in the same way we obtain for Q\ for example. 



In 



25j it is also given some direct evidence that the smoothness of the inhomogeneity profile 



is important to allow to distinguish these models from ACDM, which is consistent with our 
conclusion that the solution of the inversion problem is possible only if we allow a not smooth 
radial matter profile, while smooth models 26[ can be distinguished from ACDM, at least in 
principle, without any ambiguity 

The main point is that the luminosity distance is in general sensitive to the geometry of 
the space through which photons are propagating along null geodesies, and therefore arranging 
appropriately the geometry of a given cosmological model it is possible to reproduce a given 
Dl(z). For FLRW models this corresponds to the determination of Q\ and fl m and for LTB 
models it allows to determine the functions E(r), M(r),tb(r). 

Another observable which could be used to constraint LTB models is the redshift spherical 



shell mass mn(z) |27j . which has been recently calculated [28| for a central observer up to 
the fifth order in the redshift, and can be generalized {29 1 to the more observationally related 
E R ss{z), the redshift spherical shell energy 

Since mn(z) is not yet available observationally, we don't know if it is in agreement with 
ACDM models. It could be possible that only LTB models can fit the observed mn(z), in 
which case it would allow to distinguish between LTB and ACDM models. If on the contrary 
mn(z) is in agreement with ACDM predictions, we could still not distinguish between them, 



as shown numerically in 



30| . and analytically in this paper. In this case there would be only 



one LTB model able to exactly mimick the mathematical ACDM prediction for both mn(z) 
and Dl(z) and it would corresponds to a matter density profile not smooth at the center. 

In this paper we derive for the first time in the FRLW gauge the general formula to third 
order in redshift for the luminosity distance for a central observer in a matter dominated LTB 
space without imposing the conditions for a smooth matter distribution. 

We then calculate the low redshift expansion of mn(z) and Dl(z) for a flat ACDM model 
and show how it possible for both these observables to mimick with a matter dominated LTB 
model the effects of the cosmological constant in a flat homogeneous space, that the solution 
is unique and corresponds to a matter profile which is not smooth at the center. 

Our results are general since we don't use any special ansatz for the functions defining LTB 
models, and the arguments about the number of free parameters defining the LTB models 
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constrained by observations are independent of the value of the cosmological constant. 

We need three functions to define a LTB solution, but because of the invariance under 
general coordinate transformations, only two of them are really independent. This implies that 
two observables are in principle sufficient to solve the IP of mapping observations to a specific 
LTB model, for example the luminosity distance Dl(z) and the redshift spherical shell mass 
m(z)n(z) = mn{z). 



As observed by [30| , there has been sometime some confusion about the general type of LTB 
models which could be used to explain cosmological observations, so it is important to stress 
that without restricting the attention on models with homogeneous big bang, a void is not 
necessary to explain both Dl(z) and m(z)n(z) = mn(z) without the cosmological constant, 
but we will prove that this is only possible for LTB models that are not smooth at the center. 



II. LEMAITRE-TOLMAN-BONDI (LTB) SOLUTION 



Lemaitre-Tolman-Bondi solution can be written as 
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{R, r ) 2 dr 2 



ds 2 = -dr + 



+ R 2 dn 2 



(1) 



1 + 2E 

where R is a function of the time coordinate t and the radial coordinate r, R = R(t,r), E is 
an arbitrary function of r, E = E(r) and R, r = OR/ dr. 
Einstein's equations give 



p{t,r) 



2E(r) 2M(r) 
2M, r 



R 3 



R 2 R, 



(2) 
(3) 



with M = M(r) being an arbitrary function of r and the dot denoting the partial derivative 
with respect to t, R = dR(t,r)/dt. 
We can introduce the variables 



A (t,r) = ^-t, k{r) 



2E(r) 



Po{r) 



6M(r) 



(4) 



so that Eq. (TfJ) and the Einstein equations ()2]) and ([3]) are written in a form similar to those 
for FLRW models, 

2 dr 2 



ds 2 = -dt 2 + A 2 



1 + 



A )r r 



A J 1 — k(r)r 2 



+ r 2 dVl 2 



(5) 
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(J) 

p{t,r) 



K r ) , po( r ) 

(P0r 3 ),r 

3A 2 r 2 (Ar), r ' 



(6) 
(7) 



The solution can now be written as 



a{r],r) 



Po{r) 
6k(r) 

Po(r) 
6k(r) 



cos 



k(r) r] 



sm 



k(r) 



k(r) i] 



+ t b (r) 



(9) 



where r) = r r = / dt' /A(t', r) and A(t(?7, r), r) = a(r/, r) and t&(r) is another arbitrary function 
of r, called the bang function, which corresponds to the fact that big-bang/crunches can happen 
at different times. This inhomogeneity of the location of the singularities is one of the origins 
of the possible causal separation [U[ between the central observer and the spatially averaged 
region for models with positive an- 

In the rest of paper we will use this last set of equations. Furthermore, without loss of 
generality, we may set the function po( r ) to be a constant, po( r ) — Po — constant, corresponding 
to the choice of coordinates in which M(r) cx r 3 , and we will call this, following 34|, the FLRW 
gauge. 



III. GEODESIC EQUATIONS 



We will adopt the same method developed in [35| to find the null geodesic equations in the 
coordinates (i],t), but here instead of integrating numerically the differential equations we will 
find a local expansion of the solution around z = corresponding to the point (to, 0) = (770, 0), 
where t = t(r/o,0). We will indeed slightly change notation to emphasize the fully analytical 
r.h.s. of the equations obtained in terms of (77, t), on the contrary of previous versions of the 
light geodesic equations which require some numerical calculation of R(t, r) from the Einstein's 
equation(T2]). 

For this reason this formulation is particularly suitable for the derivation of analytical results. 
The luminosity distance for a central observer in a LTB space as a function of the redshift 
is expressed as 

D L (z) = (1 + z) 2 R (t(z), r{z)) = (1 + zfr\z)a (rj(z), r{z)) , (10) 
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where \t(z),r(z)j or \Jj](z),r(z)j is the solution of the radial geodesic equation as a function 
of the redshift. 

The past-directed radial null geodesic is given by 

dT ^ f(T( \ \ f( + \ ~ R At,r) ,,,, 
= /(T(r),r); f(t,r) = —= (11) 



dr ■ y / l + 2£(r) ' 



where T(r) is the time coordinate along the null radial geodesic as a function of the the coor- 
dinate r. 

From the implicit solution, we can write 

T(r) = t(U(r),r) (12) 
dT(r)_ = todUtf + m 
dr dr/ dr dr 

where U (r) is the rj coordinate along the null radial geodesic as a function of the the coordinate 
r. The geodesic equations for the coordinates (rj, r) can be written as, 

dr] d r t(r],r) — F(r),r) 



dz (1 + z)d v F(r),r) 

dr a(^, r ) 

dz (1 + z)d v F(r], r) 



P(v,r), (14) 
q(r],r), (15) 



F(r),r) = ===== [d r (a(r),r)r) + d v (a(ri,r)r)d r rj\ , (16) 

Jl — k(r)r 2 

where r] = U(r(z)) and F(r),r) = f(t(r),r),r). It is important to observe that the functions 
p, q, F have an explicit analytical form which can be obtained from a(r), r) and t(r), r) as shown 
below. 

The derivation of the implicit solution a(r),r) is based on the use of the conformal time 
variable 77, which by construction satisfies the relation, 

dr)(t, r) 



dt 

This means 



= a-K (17) 



n 1 1 

t(r),r) = t b (r)+ / a(r) ,r)dr) , (18) 
Jo 

dt = a(r], r)dr, + ^jH ^^'^ dr]' + t' b (r^J dr , (19) 
In order to use the analytical solution we need to find an analytical expression for F and F v . 
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This can always be done by using 



d_ 

d r t(rj,r) 



dr] 



Po k'{r) 



d r rj = —a(i],r)~ 1 d r t 



3 sin ( r)Jk(r) ) — r] ( 2 + cos ( rjyjk(r) ) \Jk{r) 



+ f h (r) 



(20) 

(21) 
(22) 



In this way the coefficients of equations (fT4l) and ( IT51) are fully analytical, which is a significant 
improvement over previous approaches. 



IV. CALCULATING D L {z) AND mn(z) 

Here we will not give the formulae in terms of 770 an d trigonometric functions, since they 
are rather complicated and not relevant to the scope of this paper, but rather introduce the 
following quantities: 



tan(^p) 2 po , x 

a = a 770, = ^ , 23 

3fc (tan(^) 2 + 1) 



3k 3 /2 (tan(^a) 2 + 1 



Ho = ^-T^f- L , (24) 



tan(«) 3 p 

, 4(tan(^ + l), (25) 



where we have used 



H = (26) 
a(t ,0) 

a(t ,0)a(t o ,0) 

q °-~ a(t Q ,0f ■ (27) 

The derivative respect to t is denoted with a dot, and is calculated using the analytical 
solution a(r],r) and the derivative respect to r\ : 

a = d t a = d v a a~ l . (28) 

In deriving the above equations we have expressed all the relevant quantities in terms of 
tan( v ^° ,?0 ) using standard trigonometric identities such as cos(x) = YTIan^Or/f)' This has the 
advantage to make formulae valid for any sign of fco by analytical continuation and to obtain 
more compact expressions, which are easier to simplify. Without such an approach the results 
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of the calculations can in fact become quite cumbersome and difficult to treat even with a 
mathematical software such as MATHEMATICA, which was actually used to derive the results 
presented in this paper. For this reason a series of simplifying routines have been developed 
based on the simplifying procedure of recursively re-expressing any trigonometric function in 
terms of tan( V ^ T?0 ), allowing to reduce to a purely algebraic operation the final simplification. 

As an extra check all the results have been independently derived using an alternative 
method based on the local expansion of the solution of the Einstein's equations in terms of the 
variables (t, r), but we will report this in a separate paper in preparation about the deceleration 
parameter in LTB models. 

Another convention we will follow in the rest of the paper will be to express everywhere 
ko, po, rjo in terms of H , a , qo by inverting equations (I23H25I) : 



k = aXtan 2 (^), (29) 
po = 3aX 2 [l + tan 2 (^)], (30) 



2 arctan ^J2% ■ 
Vo = - , „^ ■ (31) 

2 



a # tan(») 



We may actually set ao = 1 by choosing an appropriate system of units, but we will leave it 
in order to clearly show the number of independent degrees of freedom of the problem and to 
emphasize the difference with AC DM models. 

Expanding the r.h.s. of the geodesies equation we can easily integrate the corresponding 
polynomial q(z),p(z), to get r(z) and r](z). It can be shown that in order to obtain Dl{z) to 
the third order and mn(z) the fourth we need to expand r(z) to the third order and r](z) to 
the second. In this paper we will try to solve the IP using the expansion: 



k(r) = k + far + k 2 r 2 + .. (32) 
t b (r) = t b Q + t\r + t b 2 r 2 + .. (33) 



It has been proved 



36] that the energy density is smooth only if the linear term in the above 
expansion is zero, i.e. fa = t\ = 0, or the Laplacian in spherical coordinates would diverge at 
the center. This implies that including linear terms in the above expansions we are considering 
models which are not smooth at the center. For definiteness we will present here the results of 
the calculations to the second order in r, corresponding to respectively third and fourth order 
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for Dl(z) and mn(z). We will show later nevertheless that our conclusions about the existence 
and uniqueness of the solution are independent of the order at which we truncate the above 
expansion. 

It is important to observe that in general, up to the order at which we will consistently 
expand all the relevant quantities, we have eight parameters : 

p ,r}o,t b ,t b 1 ,t b 2 ,ko,k 1 ,k 2 (34) 

or equivalent ly 



-^0)90)^0)^0)^1)4)^1)^2 (35) 

where, from eq. fl23H25"]) . we have used the fact that H ,q ,a contain the same information 
of k ,r] ,p , 

There is one important question : should we regard the conformal time coordinate of the 
central observer 770 as an independent parameter? Or there are some constraint to make it 
consistent with the age of the universe, which, if we assume the inhomogeneities to be only 
local, should be approximatively the same as the one estimated in AC DM models? 



As observed in 



30j | since LTB solutions correspond to a pressureless spherically symmetric 



matter distribution, they cannot be extended to the earliest stages of cosmological history, 
implying that rjo is related to the time when this kind of inhomogeneities arise, more than 
to the age of the Universe, and since it is re-absorbed in the definition of Hq, go, a>o given in 
eq. fl23H25"]) . it is correct to consider it as an independent parameter. 

As we will show later, observable do not depends on do, so it should not be considered a 
really independent parameter. This is a consequence of the fact that H ,qQ,ao are invariant 
under changes of parameters which preserve a/^o^o- 

We get six constraints from the expansion of mn(z) and Dl(z) respectively to fourth and 
third order, so in principle without imposing the smoothness conditions t\ = k\ = we should 
be able to solve the IP of locally mapping the two observable of a LTB model to any given 
ACDM model, even taking into account the fact that, as we will show later, mn(z) and D^(z) 
don't depend on £q, and the solution should be unique. 

It is also clear that a simple preliminary argument based on counting the number of inde- 
pendent parameters implies that the IP cannot be solved if the matter distribution is smooth 
at the center, i.e. if t\ — k\ — 0, since we would have six constraints and only four truly 
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independent parameters, and this was shown more in detail in 26]. 



After re-expressing the results in terms of Ho,qo, eto and X = arctan \/2q — 1 we get: 



r?0) 
r(z) 



r l z + r 2<2 2 + r 3^ 3 



(36) 
(37) 



2fci(g + l)^2qJ=TX - (2q - l)(2a 3 H 3 (l - 2q ) 2 + 2a 2 H$(l - 2q ) 2 t b 1 + 3h 

2a 4 ,# 4 (2g - l) 3 

- J2q -lX(2a 3 H 3 k 1 (Uq 3 - 3q 2 - 1) + 



4a 7 # 7 (l-2g ) 4 

-4a 2 # 2 (2g - l)(2# &i(l - 2q )q t b 1 + A; 2 (g + 1)) + k 2 (llq + 5)) + 
+ (2g - l)(2a 6 H 6 (q + l)(2g - l) 3 + 2a 5 H 6 (2q - l) 3 (3g + l)t\ + 
+2a 4 l # 5 (2g - l) 3 (#o(2g - l)tf - 2t\) + a 3 F 3 A; 1 (26g 2 - llg - 1) + 
+2a 2 H 2 (2q - l)(2H k 1 (2q - l)t\ - 3k 2 ) + 9kj) + Qk 2 lQ 2 X 



Tl 



r 3 



a H 



v^^T (-a 3 H$ (4g 3 - 3g + 1) - 2a 2 # 3 (l - 2q Q ) 2 q Q t\ - 5hq + h) + 6hq 2 X 

2a^ 4 (2g - If' 2 
6(2g - l)g 2 X(8a 2 ^ 3 A; 1 (2g - l)(a g + (2q ~ 1)4) + 



8alHZ(2q - If/ 2 

+4a 2 H 2 k 2 a - 2q ) + 17k 2 ,) + (2q - lf 2 (4a e H e (2q - l) 3 (g 2 + 1) + 16a^ 6 g 2 (2g - 1)H\ + 
+8a 4 l # 5 g (2g - l) 3 (H (2q - l)tf - t\) + 4a 3 # 3 £; 1 (i2g 3 - 2g 2 - 4g + 1) + 
+4a 2 H 2 (2q - l)(6H k iqo (2q ~ ^4 ~ ^ 2 q + k 2 ) + A; 2 (43g - 5)) + 



+72k\ J2q - lq 3 X' 



We can then calculate the luminosity distance : 



Dl(z) 
D 2 



(1 + z) 2 r{z)a{r]{z), r(z)) = D lZ + D 2 z 2 + D 3 z 3 + .. 
1 

1 

2a 3 M2q - If/ 2 
+Qk iqo {q + 1)X 



(3? 



2q - l(-a 3 # 3 (l - 2g ) 2 (g - 1) - 2a 2 # 3 (4g 3 - 3g + l)t\ - 9fc l9o ) + 



4Q 6 g 7 (2go _ 1)11/2 [_ - 3(2g - l)g X(4a 2 F 3 A; 1 (l - 2q ) 2 q (2a q + Aq Q t\ + t\) + 
-Aa 2 H 2 k 2 (Aq 3 - 3g + 1) + A; 2 (50g 2 + 7q - 7)) + (2q - l) 3 / 2 (2a^ 6 (l - 2g ) 4 (g - l)?o + 
+8a^ 6 (l - 2q ) A q 2 t b 1 + 4a 3 # 3 £; 1 (l - 2g ) 2 go(5g - 1) + 

+2a 2 # 2 (l - 2q ) 2 (H k 1 (20q 2 + q - l)t\ - 9k 2 q ) + 2# 5 (a - 2a q f(H q (Aq + l)(tf + 
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-2(q + l)t b 2 ) + 9A^g (Hgo - 4)) + 18^2% - l(4g + 1)% 3 ^ 2 

The effects of inhomogeneities are showing already from the second order, because we are 
not setting the smoothness conditions t\ — k\ — 0. 

From the definition of mn(z) and the equation for the energy density we can write 



Anmn(z)dz = pd 3 V 
from which by using dr = (dr/dz)dz we get 



8nM' 



-.dr 



mn(z) 



2M'(r(z)) dr(z) 



yjl — k(r)r- 



p r{z) 



(39) 



dr(z) 



'I - k(r(z))r(z) 2 dz yj\ - k(r(z))r(z) 2 dz 
where in the last equation we have used the FLRW gauge condition M(r) 
allows to calculate mn(z) directly from r(z). 
We finally get: 



(40) 

Por 3 /6, which 



mn{z) 
mn 2 



mn 2 z 2 + mn^z 3 + mn^z 4 
6g 
H 



a 3 M2q - 1)V2 
■\2kx<ftX- hq 2 sm(4X)) 
2 



(41) 
(42) 

3q (2J2q - l(a 3 # 3 (4g 3 - 3g + 1) + 2a 2 # 3 (l - 2 9o ) V? + *hqo) + 



(43) 



777,714 



32a[|#o 7 (2go - l) 1 



3g (120g 2 X( v /2g - l(Wi(l - 2g ) 2 (3g + 1) + 
-4a 2 (# - 2H q ) 2 (2H k 1 (l - ?>q Q )t\ + k 2 ) + k 2 (50q - 17)) - 2k\ql sin(4X)) + 
-4(2g - l)(2a^ 6 (l - 2g ) 4 (15g 2 + 14g + 13) + 40a^ 6 (l - 2g ) 4 g (3go + 1)*? + 



+20a 3 # 3 A; 1 (l - 2q Q f(\\q 2 + 6g - 2) + 20a 2 # 2 (l - 2q ) 2 (2H k iqo (llq - A)t\ + 



-5k 2 q + k 2 ) + 40H*q (a - 2a g ) 4 (#o(3g - l)(t?) 
-UA0k 2 (3q - l)q 3 X 2 ) 



4) + 5fc 2 (136g 2 -73g + 7)) + 



On the contrary to the smooth case, it can be seen that the effects of inhomogeneities show 



already from the third order, and while also the fifth order has been calculated [28j , we don't 
report it here since we don't need it. 

A final important point is that the two observables don't depend on t\, which implies that 
we are left with six parameters to determine in order to solve the IP, while we will have six 
equations for the coefficient of the expansions. 
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V. CALCULATING D L (z) AND mn(z) FOR KCDM MODELS. 



The metric of a kCDM model is the FLRW metric, a special case of LTB solution, where : 

Po(r) = Po (44) 

k(r) = (45) 

t b (r) = (46) 

a(t,r) = a(t) (47) 

In this section we will calculate independently the expansion of the luminosity distance and 
the redshift spherical shell mass for the case of a flat KCDM. 

We can also use these formulas to check the results given in the previous section, since in 
absence of cosmological constant and inhomogeneities they should coincide. 

Here we will denote all relevant physical quantities with an upper script A to avoid confusion 
and clearly distinguish them from the quantities introduced previously for LTB models. 

One of the Einstein equation can be expressed as: 



,A\3 



h\ z ) = h a ^ n m +n A = H A ^/n m (i + z) 3 + n A (48) 

We can then calculate the luminosity distance using the following relation, which is only 



valid assuming flatness: 



D A (z) = (1 + z) I = D A z + D A z 2 + D A z 3 + ... (49) 



From which we can get: 



D? = 77a (50) 



A An A + n 



A = ^ , u uM 



= ^ (52) 
3 8# A K ' 

It is convenient to re-express the above coefficient in terms of the two observables H A , q A 

D i =W (53) 

,A\2 , „A 



,A _ % A ) 2 + q A - 2 
QH A 



D A = (55) 
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where we have used the following relations 

n L + n M = i (56) 
n M = (57) 

It should be underlined here that Hq, a A appearing in this formulas are not the same as 
the ones defined in the previous section for LTB models. 

For calculating mn(z) we first need r(z), which can be obtained from the radial null geodesic 
equation which in this case takes the simplified form 

dr A {z) 1 1 



dz ' (l + z )a A a A H A 

and after integration we get 



(58) 



1 r z dz' 

U = hkw&= (59) 



a A Jo H A (z') 

(g A + l)z 2 (3(^ + 4^ + 1)^ 



a A H A 2a A H A Qa A H 



(60) 



We can now calculate mn A (z) 

p = 3(a A ) 3 Q M (H A ) 2 (61) 

dr A (z) 

mn A (z) = mn A z 2 + mn A z 3 + mn A z A = p r A (z) 2 v = (62) 
2{q A + l)z 2 4{q A + l) 2 z 3 5{q A + l) 2 {9q A + 5)z 



+ ° W0 ' 1 ^7°'^ + O (, 5 ) (63) 



H A H A QH A 

We can check the consistency between these formulae and the one derived in the case of 
LTB without cosmological constant by setting: 

k x = k 2 = t\ = t\ = (64) 
g = 1/2 (65) 

which corresponds to the case in which Q M = 1. 

It is important to mention again that a A ,q^, Hq defined in this section are in general different 
from the ones defined in the previous section for LTB models. 
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VI. LOCAL SOLUTION OF THE INVERSION PROBLEM 



In this section we will denote with an upper script A all the relevant quantities referred to 
a AC DM model, including the coefficients of the redshift expansion for Dl(z) and mn(z). 

In order to solve he inversion problem to the second order in r we need to solve the following 
system of six equations equations: 

Di=Di 1 < i < 3 (66) 
mrij = mrtj 2 < j < 4 (67) 

Including higher orders in the expansion would not change our conclusions since if we include 
a new order for k n ,t b n we would than have two new parameters to solve for and two new 
conditions coming from the matching of the coefficient {D n+ i, mn n+ 2] in the redshift expansion. 

We can first solve the system of two equations 

= D 1 (68) 
mn* = mn2 (69) 

which gives the simple solution: 



H = H£ (70) 
* - ^ (71) 

Since q$ > —1, the above equation implies q > 0, so that the positivity condition implied 
by eq. (!25j) is respected, and the inversion problem is possible. We are left with four parame- 
ters {ki, &2, t\, t\} to determine, and since we have only four more independent equations, the 
solution should be unique, but we will report it in detail in a separate work, while here we are 
interested in proving the existence of such a solution. 

Another important consequence of our calculation is that it is impossible to solve the in- 
version problem for both mn(z) and Dl(z) setting one the two functions k(r) or t b (r) to zero, 
even allowing not smooth matter profiles, since we would then have less parameters to solve 
for than equations. 

On the other side as long as we are trying to solve the inversion problem for a single observ- 
able, we would always have more free parameters than equations, even imposing the smoothness 
conditions, meaning that in this kind of IP there is not a unique local solution. We could use 
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this freedom for example to smooth the density profile. We will further investigate this a in a 
future work, where we will study the local solutions to the IP for D t,(z) and mn(z) separately 
and for both the observables at the same time. In the smooth case 26] there is no solution of 
the IP for both Diz) and mn(z), since there are not enough free parameters to solve the IP 
equations. In the smooth single observable IP for Di(z) we obtain instead go = 1o, which im- 
plies the well known fact that negative g^ cannot be reproduced by LTB models with smooth 
matter distribution |36|], since go is always positive. This single observable smooth case IP 
solution is completely different from the not smooth multiple observable case in eq. (]71j) since 
they correspond to different types of IP, i.e. they satisfy different set of equations. For this 
reason the single observable smooth solution cannot be obtained from the multiple observable 
solutions in the k\ — t\ = limit because they are different sets of solutions, and in fact the 
IP has no solution in the multiple observable smooth case. 

It should be also mentioned that since according to our analysis we are left with no free 
parameter after solving the IP for both mn(z) and Dl{z), it is highly probable that additional 
independent observables should make possible for a central observer to distinguish without any 
ambiguity between matter dominated LTB and AC DM models. 



VII. DO WE NEED A VOID TO EXPLAIN LUMINOSITY DISTANCE OBSERVA- 
TIONS? 

So far we have focused on the possibility to solve the IP for both mn(z) and Dl(z) in 
presence of a centrally not smooth spherically symmetric matter distribution, showing how the 
on the contrary to the smooth case a solution should exist, as far as Hq and go are concerned, 
giving they explicit relation with qfr and Hq\ 

ttLTB ttA „LTB 1 + 1o / 70 \ 

H o ~ H o > % = — 3 — • I 72 ) 

It actually turns out that the solution is unique, since after introducing the two following 
dimensionless parameters : 

K x = (a tfor 3 fci , K 2 = (a H )- i k 2 (73) 
Tl = a^t\ , T* = av l (a H o yHl (74) 

we can eliminate in the formulae for the observables, which means that, as expected, they 
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do not depend on do- For example we can get : 

6ir lg o( go + 1)V2^TX - (2g - 1) (gpg^i - 6T, 1 + 5) + ^(ST, 1 + 4) - 8g 2 + 2T, 1 - 1) 
2 2H (2q - 1)3 

(75) 

This can be seen as a consequence of the fact that i]q and k always appear in the form y/korjo 
in the definition of a ,q ,H , implying that they do not really correspond to two independent 
degrees of freedom of the problem, since physical quantities are invariant under changes of 
parameters which preserve a/^o^o- Since £° never appears in the expansion of the observables, 
we are left with only six independent parameters and six independent constraints, implying 
that there is at most only on independent solution to the IP 

A detailed analysis of the solution is left to a future work, while here we will consider the 
consequences of our results on fitting the luminosity distance data. Let us start from considering 
the constraints coming from the first two terms: 

D = d£ , D 1 = D^ (76) 

which imply that out of the four independent parameters H Q ,q , Ki,Ti, two of are left free, so 
the IP is degenerate, even considering a homogeneous bang function, i.e. setting T\. This is 
confirmed by the fact that several different solutions have been prosed in the literature with 
tb{f) — 0. Going to a higher order would increase the degeneracy of the solution, since we 
would then only have one extra constraint coming from D n and two additional parameters to 
determine, K n ,T n . 

To give a more definite answer to the question posed in the title of this section, we should 
actually look for the best fit parameters and then calculate the corresponding density profile, 
but this would go beyond the scope of the present paper which is to prove that a centrally not 
smooth matter distribution can mimick the cosmological constant for both mn(z) and Drz). 

We can nevertheless observe, based on the recent numerical work by [30j and our analysis 
performed so far, that in presence of a not homogeneous bang function there exist solutions 
of the IP for both mn(z) and Dl(z) which correspond to overdense regions. This solutions 
are obviously valid even in the case in which we consider only Dl(z), implying that as far 
as D L (z) is considered, both overdense or underdense regions can fit the data if we include 
inhomogeneous bang functions. 

In the homogeneous bang case t&(r) = 0, only solutions corresponding to a central void are 
known, and we leave to a future work a detailed study of this case based on the low redshift 
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expansion we have derived, in order to determine if this is the only type of possible models, but 
based on physical arguments about the dynamics of void expansion we expect that the dimming 
of supernovae light can only be explained by an underdense region for the homogeneous bang 
case. 

One last important comment is that so far we have focused on reproducing the coefficient 
of the redshift expansion of the observables corresponding to the best fit flat ACDM models, 
but this kind of approach is not completely rigorous since these models depend only on the 
two independent parameters Hq,Qm, while LTB models have a higher number of parameters, 
implying that they could actually provide an even better fit of experimental data. 

VIII. CONCLUSION 

We have calculated for a central observer in a LTB space without cosmological constant 
the low redshift expansion for the luminosity distance and the redshift spherical shell mass 
respectively to third and fourth order. We have then derived the same observables for ACDM 
models and shown how it is in principle possible to locally mimick a ACDM model with a LTB 
matter dominated universe if the conditions for a centrally smooth matter distribution are not 
imposed, i.e. if we don't set k\ — t\ — 0, and the solution is unique. 

Another important consequence of our calculation is that it is impossible to solve the IP for 
both mn(z) and D L (z) setting to zero one the two functions k(r) or t b (r), even allowing non 
smooth matter density profiles, since we would then have more independent equations to solve 
than parameters to determine. 

For the case of the inversion problem for a single observable we have an even wider space of 
solutions and we could try to use this freedom to impose the conditions for a smooth matter 
profile. It is in fact possible to find a smooth matter distribution able to reproduce the best fit 
ACDM redshift expansion for the luminosity distance if the bang function is not homogeneous, 
and we will investigate this in detail in a separate work. 

It should be also mentioned that since according to our analysis we are left with no more 
free parameters after solving the IP for mn(z) and Dl{z), it is highly probable that additional 
independent observables would make possible for a central observer to distinguish without any 
ambiguity between matter dominated LTB and ACDM models. 

We will further investigate this in a future work, where we will study in detail the local 
solutions to the IP for D L (z) and mn{z) separately, for both the observables at the same 
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time, and for other relevant observables, determining when overdense or underdense regions 
are possible solutions. 

One last important comment is that we have focused on reproducing the coefficient of the 
redshift expansion of the observables corresponding to the best fit flat ACDM models, but this 
kind of approach is not completely rigorous since these models depend only on the two inde- 
pendent parameters H , Qm, while LTB models have a higher number of parameters, implying 
that they could actually provide an even better fit of experimental data. From an observational 
point of view it would in fact be more important to confront directly with actual experimental 
data rather than reproduce the best fit theoretical ACDM model, so our conclusions should be 
considered keeping this limitation in mind, and we leave to a future work such a direct analysis 
of experimental data. 
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